Effect of degree correlations on the loop structure of scale-free networks 
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In this paper we study the impact of degree correlations in the subgraph statistics of scale- 
free networks. In particular we consider loops, simple cases of network subgraphs which encode 
the redundancy of the paths passing through every two nodes of the network. We provide an 
understanding of the scaling of the clustering coefficient in modular networks in terms of the maximal 
eigenvector of the average adjacency matrix of the ensemble. Furthermore we show that correlations 
affect in a relevant way the average number of Hamiltonian paths in a three-core of real world 
networks. We prove our results in the two-vertex correlated hidden variable ensemble and we check 
the results with exact counting of small loops in real graphs. 
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I. INTRODUCTION 



The dynamics and the function of many complex sys- 
tems strongly affect their network structure 00, E| • 
In fact both large-scale properties (like the scale-free de- 
gree distribution la) and local properties (like recurrence 
of small motifs |g, Q) must be selected for widespread 
robustness requirements and specific preferential uses in 
real graphs. A large number of different networks 0,0,01, 
from the Internet to the protein interaction networks in 
a cell, share a scale-free degree distribution P{k) ~ k^ 1 
with 7 < 3 and a high clustering coefficient respect to 
random Erdos-Renyi graphs Q. The scale-free degree 
distribution of a network affects the statistics of sub- 
graphs present in it showing that large-scale properties 
and local properties of scale-free networks are strongly 
related to each other. Special examples of subgraphs in 
networks are loops [(J, Il3 > paths that pass through each 
node in the loop only once. In random scale- free networks 
there are many small size loops compared to random 
graphs and there can be a lack of Hamilton cycles (loops 
of length L = N) due to the fact that most of the large 
paths need to pass through hubs 10 . Along with other 
properties, many real scale-free networks also have degree 
correlations Degree correlations in real networks in- 
dicate that links are not randomly wired and that the 
probability that two nodes of degree fcj and kj are linked 
deviates from the expected value r,-j = kikj/((k)N). 
Consequently, correlated networks have at least one of 
the three following features: i) a k dependent average 
connectivity k nn (k) of the first neighbors of a node with 
degree k 0,0; ii) a non trivial dependence on the con- 
nectivity of the clustering coefficient C = C{k) of nodes 
of degree k [l4|; Hi) a cutoff that is larger than the struc- 
tural cutoff K ~ {k)N. In particular many real scale- 
free networks show a power-law dependence on k both for 
k nn (k) and for C(k), i.e. k nn (k) ~ k a and C{k) ~ k~ s . 
Correlations do affect the subgraph statistics as shown in 
the Internet 0] and in calculations based on the scaling 
of the clustering coefficient [7|, [16J . Every network can be 
represented in terms of its adjacency matrix ((a)) of ele- 
ments Q>i,j = 1,0 depending if there is a link between node 



i and node j. From a formal point of view an ensemble of 
networks is given when a probability V(a) is assigned to 
each adjacency matrix ((a)) of Nx N elements. In an un- 
corrected and undirected network ensemble with given 
degree sequence {k{\ all the links are independent. Con- 
sequently all the matrix elements a^j with i < j of are in- 
dependent and their average value in the ensemble can be 
written as (cnj) = r^j — • A two- vertex correlated 
network is a network in which still the matrix elements 
cii j with i < j are independent but (a,^-) = r,-j 7^ ^k)N • 
Networks with higher order correlations instead would 
have non independent matrix elements which will favor 
some specific motifs in the network. In this paper we are 
going to provide an analytic calculation of the number of 
loops in two- vertex correlated scale- free networks. In the 
light of our results we are able to interpret the scaling of 
the clustering coefficient C(k) in terms of the scaling of 
the maximal eigenvector (the eigenvector associated with 
the maximal eigenvalue) of the average adjacency matrix 
of the network ensemble. Moreover we show that the 
maximal eigenvalue and the corresponding eigenvector 
not only determine the number of triangles in the two- 
vertex correlated network ensemble, but they also fix the 
number of small loops of length 3 <C L <C N . Finally we 
are able to give a sufficient condition for the absence of 
Hamilton cycles in two- vertex correlated networks. This 
allows us to study a set of real graphs (the Internet at 
the Autonomous System -AS- Level and protein-protein 
interaction networks-DIP) an d show that assuming 
they are specific instances of two-vertex correlated net- 
work ensembles, one can exclude the presence of Hamilto- 
nian cycles in the three-core of these graphs. Our findings 
are in agreement for the Internet with what was found in 
Ref. |l8j where a Belief-Propagation algorithm was ap- 
plied to the measurement of the number of loops in real 
graphs. The absence of Hamiltonian cycles in a three- 
core of a network is an unexpected result since regular 
random graphs with connectivity c > 3 are Hamiltonian 
[8|,[l!j. We note here that the average number {Ml) °f 
loops of size I in a two-vertex correlated network can 
possibly be dominated by a very large number of loops 
occurring in very rare networks 0]. Nevertheless, pre- 
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liminary results indicate that in uncorrelated scale-free 
networks with 7 < 3 the ratio nj^h is bounded at least 
for small loops and for Hamiltonian cycles. We expect 
that similar arguments could also be extended to scale- 
free correlated networks. 

The paper is organized as follows: in Sec. II we give an 
intuition of the results found for small loops in two- vertex 
correlated network ensemble by considering the problem 
of exact counting of loops in generic networks; in Sec. Ill 
we introduce the hidden variable ensemble and we cal- 
culate the number of small loops and Hamiltonian cycles 
in two vertices correlated hidden variable ensembles; in 
Sec. IV we compare the results with real networks; and 
finally we give the conclusions in Sec.V. 



where the factor 2L accounts for the multiplicity of nodes 
a single loop pass through and the two possible directions 
of each loop. Thus we found by intuitive arguments that 
the total number of small loops of size L of scale- free net- 
works will scale like Ag while the number of small loops 
passing through a node is proportional to the square of 
the maximal eigenvector associated with Ao. These argu- 
ments apply for the exact counting of small loops in real 
networks. In a random graph ensemble the adjacency 
matrix is a random variable which has average values of 
the elements (a,j) = r^j and we need to evaluate the 
average number of loops {Ml) instead of Ml The results 
we will prove in the following sections are an extension 
of the expressions @J and J2J to two-vertex correlated 
hidden variable network ensemble. 



II. COUNTING SMALL LOOPS IN REAL 
NETWORKS 

In this section we would like to provide some intu- 
itive arguments to show that in scale-free networks the 
maximal eigenvalue of the adjacency matrix and the cor- 
responding eigenvector are responsible for the number of 
small loops present in it. The adjacency matrix ((a)) of 
a simple network of size N is the N x N matrix of el- 
ements ctjj = 1,0 indicating the existence {a,ij — 1) or 
not (cii.j = ) of a link between node i and node j. The 
total number of closed paths of length L passing though 
a node i is given by the matrix element (a L )i^. The loops 

M L of size L passing through a node i are given by 



Ml = (a L )i,i — (corrections) 



(1) 



where these corrections account for closed paths which in- 
tersect themselves at least once and which must be sub- 
tracted from the term (a L )i } i in order to consider only 
loops. If by A„ we indicate the eigenvalues and by u n 
the eigenvectors of the adjacency matrix ((a)) we find |9( 



A/f 



(n) (n) 



£>(u^ 4 ) 



(2) 



For small L, the correction terms can be neglected if the 
spectrum of the graphs {A} contains one large eigenvalue 
Ao = Ao and if the associated normalized eigenvectors 

(n) 

satisfy < u\ <C l,Vi, as is the case in most scale- free 
networks. If these conditions are satisfied the sum over n 
in is dominated by the term n = and consequently 
the number of loops of length L passing through the node 
i is given by 



Ml 



(i) 



Aq UJ'II 



o„,o 



(3) 



while the total number of loops of size L is given by 



(4) 



III. AVERAGE NUMBER OF LOOPS IN 
CORRELATED HIDDEN VARIABLE ENSEMBLE 

To model a general two-vertex correlated network in 
the following we will consider networks that are gener- 
ated within the hidden variable model [20L |2*H . The pre- 
scription of Ref. [2(| to generate a class of scale-free 
networks with exponent 7 is the following: 1) assign to 
each node i of the graph a hidden continuous variable 
qi distributed according to a p{q) distribution. Then 2) 
each pair of nodes with hidden variables q, q' are linked 
with probability r{q,q'). When the hidden variable dis- 
tribution is scale-free p(q) = poq 1 for q £ [m, Q] and 
r { < lil') ~ Qi' I ((°) N) , we obtain a random uncorrelated 
scale-free network. In this specific case a structural cut- 
off is needed to keep the linking probability smaller than 
one, i.e. Q 2 /{{q)N) < 1. This cutoff scales differently 
with the system size N depending on the value of 7: 
Q ~ AtV(t-i) for 7 > 3, Q - N 1 / 2 for 7 e (2,3) and 
Q ~ jV 1 / 7 for 7 S (1, 2). On the contrary, to generate a 
correlated scale-free network with natural cutoff N 1 / (t -1 ) 
and 7 > 2 in the literature different ansatz have been pro- 
posed [20I l2l| . In order to present general results on the 
average number of loops in the hidden variable ensemble 
for any type of linking probabilities r(g, q') we consider 
an ordered set of distinct nodes . . . ,i n , . . . ,il\- With 
each such kind of set it is possible to associate a loop in 
the network in which subsequent nodes are linked with 
each other. For each choice of the nodes . . . ,%l}, 
with hidden variables {qi t , . . . , q lL } the probability that 
they are connected in a loop is 



r{qn , qi 2 )r(q i2 , q is ) ■ ■ ■ r(q ih , q h ) = r(q in , q it 



od(n + l,L) I 



(5) 

and for each loop of the network there are 2L ordered 
sets . . . ,%l} which describe it corresponding to cyclic 
permutations of the indices and to their order inversion. 
The average number of loops of size L in the graph is 
given by the number of ways we can choose an ordered 
set of L nodes . . . , il} multiplied by the probability 
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that these nodes are connected in all distinguishable or- 
derings and divided by 2L. In order to proceed with the 
calculation, we lump together nodes with hidden vari- 
able qi € [q, q + Aq), where Aq is a small interval of 
q. In each interval of q there are N q ~ NP(q)Aq nodes 
of the network. For each choice of the L nodes, let n q 
with Y] n q = L be the number of nodes in the loop 
with qi n G [q,q + Aq) . The ways we can choose them 
within the N q nodes of the network, is given by the bino- 
mial N q \/\n q \(N q — n q )\]. Moreover let n q>q > indicate the 
nodes of a hidden variable q of the loop linked with a sub- 
sequent node of hidden variable q' in the fixed direction 



of the loop. We note that the way to choose {n qtq r} is 
given by the multinomial n q \j FJ , n q,q' ! and that the par- 
tition {n q>q r} must satisfy the conditions J2 q , n q-q' = n q 
and n q q * — n q i . Finally the number of ways in which 
one can permute the L nodes keeping n q . q i constant is 
given by Y[ q n q-- Considering all this and that the prob- 
ability Eq. JSJ that the selected nodes are connected in 
the chosen order can be written as n g9 /r(g, q') n i,i' , we 
get the following expression for the average number of 
loops (Ml) of size L. 



N q \ 



2L ^ xx„ !(JV, 



n 



n„ 



l[r( q ,q') n ™' 



(6) 



9,9' 



where the sums Yl{ n }'X){n /} are extended over all 
{n q } and such that J2 q n q = L, J2 q > n q,q' = n q 

and n q , q i — n q i and the factor 2L accounts for the 
multiplicity in which we count each loop. Introducing 



the constraints J2 q n q — L and ^2 q n q q i — n q > by ex- 
plicit delta functions, using their integral representation 
we find 



w = tl dx E^n n , (jV g i n) , nve— / v Xq ne«.-. e n^n^^'e-^ 

^J-oo {nq} q n q .{IV q n q ). g g {n q . q '} 9' 9)9 qq i 



r 



where the T>x q indicates Y[ q dx q , and the sum over {n q q r} Consequently, performing the multinomial summations 
is performed over all {n q ^} such that V ' "-9,9' = r V over {"9,9'} we S et tne following expressions: 



TO = 



2L 
1 

2L 



dx e Lc 



En 



N l 



{nj 9 

dx e - x: n 



n g !(AT 9 -n 9 )! 
n,!(7V g -n q )\ 



\J°° Dx 9 I]e v « ^r( M ')e-V 



(7) 



r 



with 



In ( ^r( g , 5 ')e 



(8) 



Notice that in Eq. (JJJ) one can safely take the limit Aq — > 
and that the average over the P{q) distribution is taken 
assuming that we focus on the limit N — > 00. In what 



follows, we will evaluate Eq. J7J in different ranges of 
L in the limit N — > 00. Assuming L > 1 we evaluate 
the integral over the variables {x q } by the saddle point 
equation finding 



r(q',q) 



Eq" r(q>,q»)e- x «" 



(9) 
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If we indicate by S q < the sum S q < = Y^q r {l' il) e Xi i we 
can cast the solution in the following form, 



(10) 



'E,' n ql r(q,q')/S' q ' 
This provides the self-consistent equation for {S q } 



(ii) 



It is easy to check that {S q } satisfying the equation 



is a solution of the Eq. 1|11|) . Inserting a delta func- 
tion S (^S q — J2 q ' n q' r { c l^ ( l l ) I Sq^j and assuming that the 
Jacobian of this transformation is 1, i.e. assuming 



S q >r(q,q') 



(13) 



and using the Stirling approximation for the factorial n q , 
the integrals over x q calculated at the saddle point take 
the values Sq nq e~ nqlnt - riq " ,+nq and the average number of 
loops of size L can be expressed as the following: 



J 



/oo r r ' 

dx e L ^ jVSq JVwqY^Xl 



■(e Sq ) exp 



n q q 



n q \{N q - n q )\ 



^Wq ( Sq - U q' J ' ) \ U ) 



Finally, performing the summation over {n q } we get 



(Ml 



^- J dx e L( - x ~V YlJvSg J Vw q exp IjV ^ln(l + e - x S q 2 exp 



-Nj2w q ,r(q,q')/S q 



N^VqSq 



where {) q indicates the average over the distribution of evaluate the saddle point equations, finding 
the hidden variables N„. In the limit JV > L > 1 we 



Sq = N 



■ r(q, q')Sq,e-* exp (- r(q', q")w q n/S q ) 



1 + S$e-* exp (- Y, q " r(q',q")w q „/S ql 
w q = -P{q) 



(2S q + j: q ,r(q,q')wq / )e~ x exp 


[-J2q' K<?> q')Wq'/Sq^ 


1 + S^e~ x exp 


f-Eg' r(q,q')w q , / S q 





S\e x exp r(q,q')w q , /S q 

1 + S\&- X exp ( - J^q' r (q, q')Wq'/S q 



(15) 



with I — L/N. In order to solve these saddle point equa- With this assumption we can rewrite the saddle point 
tions we make the ansatz equations (|15fl as 

N^2r(q,q')w q , =vS q . (16) I ' r {q,q')S ql e- x - v \ 



S q = N- 



1 + Sle- X - V 



w q = -{2 + V )P{q)- L% ' 



1+S*t 
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FIG. 1: Normalized number of triangles (empty triangles), 
quadrilaterals (filled squares) and pentagons (filled circles), 
passing through nodes of connectivity k. Data are shown 
for the Internet at the Autonomous System Level in Novem- 
ber 1997 (bottom), in the s. cerevisiae protein interaction 
database (center) and in the h. pylori protein interaction 
(top) [rHl . The solid lines indicate the predictions Mh{k) oc 5*| 
where Sk is the maximal eigenvector of the correlation matrix 
N k ir(k,k'). The data are shifted to improve the readability 
of the graph. 



which can be solved and define the value of v, v = — 1. 



Network 


2{ln(S* fc ))- 


2<InOS*)>- 




{hi(p))/N 


(Hp r ))/n 


AS 


11-97 


-4.73 


2.98 




4-98 


-5.22 


3.06 




7-98 


-5.35 


3.03 




10-98 


-5.56 


3.01 




1-99 


-5.74 


3.07 




4-99 


-6.06 


3.09 




7-99 


-6.28 


3.07 




10-99 


-6.55 


3.06 




1-00 


-6.75 


3.07 




4-00 


-7.20 


3.01 




7-00 


-7.30 


3.03 




10-00 


-7.46 


3.01 




1-01 


-7.428 


3.01 




3-01 


-7.73 


3.00 


DIP 


s. cerevisiase 


-6.46 


3.99 




h. pylori 


-4.5 


3.8 




c. elegans 


-0.66 


2.89 



TABLE I: In the table we report the value of 2(ln(S 9 )) - 
(\og(p))/N with S q satisfying Eq. il'Al assuming as the max- 
imum likelihood assumption that all the qi = ki on the 
nodes of the three-core of the Internet graphs and on the 
graphs of protein interactions \lt\. We compare the value 
of 2(ln(S q )) — {\n(p))/N calculated with the two-vertex cor- 
relation assumption on real graphs or simply assuming the 
minimal assumption r(q = k, q' = k') = 1 — e~ kk /l ~ k)N 
2(\n(S q )) - (\n(p R ))/N. We observe that real correlations 
are essential to predict the absence of Hamiltonian cycles in 
these graphs. 



A. The uncorrelated case 

In the uncorrelated case, when r(q, q') — ^ N we found 

S q = 1\J~t^ which satisfies hypothesis 1|13|) . The results 
found in this limit are the same as the ones found in |10| . 

B. Small loops 

The limit of small loop size is the limit x ^> 1. In this 
limit the saddle point equations (|17[) reduce to 

s q = Y J N i<^i') s <i'^ x+1 

q' 

w q = -P{q)r{q : q')S q e- x+1 
i = (S 2 q ) q e- X+1 - (18) 

The first equation indicates that S q is the eigenvector of 
the average adjacency matrix N q /r(q,q') with eigenvalue 
A = e x ; the second equation defines the linear relation 
between w q and S q , and the third equation fixes the nor- 
malization constant for the eigenvector S q . In this limit 
the average number of loops of size L is given by 

N L ~ ^-{K) L (19) 



where A is the maximal eigenvalue of the average adja- 
cency matrix N q ir(q,q'), with the results valid until 

*«lr (20) 

where S q is the eigenvector of matrix NP(q')r(q, q 1 ) cor- 
responding to the maximal eigenvalue A 3> max5^. We 
observe that the vector Si = S qi with i — 1, . . . , N is the 
eigenvector of the matrix ry = r(qi, qj). In other words 
{Si} is the eigenvector of the average adjacency matrix 
of the networks in the ensemble (ay) = ry. This result 
provides the extension of the arguments of Sec. I , Eq. 
(QJ to the two-vertex correlated network ensemble. 

C. Small loops passing though a given node 

From expression l|15fl one can also derive the number 
of small loops passing through a given node. One can 
easily show that 

~ ^Afe 1 (21) 

where S q is the maximal eigenvector of the matrix 
N q ir(q,q') normalized in such a way that (S q ) = IK. 
This provides the extension of the arguments of section 
I Eq. @ to a two-vertex correlated network ensemble. 
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D. Hamiltonian cycles 

The Hamiltonian cycles of a graph are loops of size 
L = N. From Eq. we find that when L = N 

the expected number of Hamiltonian cycles goes to zero 
exponentially with N if 

2(ln(5,)) < 1 (22) 

with S q satisfying 

S q =J2r(q,q')N q ,/S' q . (23) 
1' 

Consequently, in the thermodynamic limit, since 

P(M L > 0) < {M L ) , (24) 

(12211 if a sufficient condition for excluding the presence of 
Hamiltonian cycles in the network. 

IV. COMPARISON WITH REAL DATA 

To test our calculation on real graphs and forecast 
some results regarding the existence or not of Hamilto- 
nian cycles we have to assume that the real networks 
under study are a particular instance of a two-vertex 
correlated hidden variable network ensemble. Since the 
average connectivity k(q) of a node depends only on its 
hidden variable the minimal assumption one can make 
to fit real networks with the hidden variable model is 
that the average degree is a one-to-one map to the hid- 
den variable q. In this assumption maximum likelihood 
considerations force us to assume that each real graph is 
a random realization of a two- vertex correlated networks 
with q { = ki and r(q = k,q = k') = ^ k)N k N k ^ Nkl where 
Nf.f.i are the total number of links between nodes of de- 
gree k and k' and Nj~ and Nk' are the numbers of nodes 
with degree k and k' . 

This results give a very interesting interpretation of 
the dependence of the clustering coefficient on the con- 
nectivity k, i.e. C(k) ~ fc ( fc 1 _ 1 ^ A^ k ySl where Sk is the 
eigenvector associated with the maximal eigenvalue A of 
the matrix Nq>r(q,q'), in agreement with the intuitive 
arguments of Sec. I. Moreover, one can predict if in the 
three-core of the considered graph there are no Hamilto- 
nian cycles by evaluating if the condition (1221) is satisfied, 



[1] R. Albert and A.-L. Barabeasi, Rev. Mod. Phys. 74, 47 
(2002). 

[2] S. N. Dorogovtsev and J. F. F. Mendes, Evolution of 
Networks (Oxford University Press, Oxford, 2003). 

[3] M. E. J. Newman, SI AM Review 45, 167 (2003). 

[4] R. Pastor-Satorras and A. Vespignani, Evolution 
and Structure of the Internet (Cambridge University 



i.e. if 

2{\n(S q ) - ^(ln(p)» < 1 with S q = '£r(q,q')N( q ')/S q . 

where the \n(p)/N = (ln[l — (1 + q + q 2 /2)e~ q ]) corrects 
for the probability that the network in the ensemble con- 
tains nodes of connectivity k < 3 as described in [ic| . 
In particular one can compare the value of 2(ln5 g ) cal- 
culated by solving |2"5|> with r(q = k, q 1 = k') extracted 

from the data [r(q = k,q' = k') = ( k ) N % k ^ Nk/ ] with the 
value of 2(hiSq) in the simplest example of a correlated 
ensemble, i.e. the static network ensemble [2^ defined 
with r(q = k, q' = k') = 1 - exp[--7#^]. We found 
as reported in Table [I] that the real degree correlations 
are such that the presence of Hamiltonian cycles in the 
three-core of the network is very unlikely. 



V. CONCLUSIONS 

In conclusion we have evaluated the number of loops 
of any size in two- vertex correlated networks. The re- 
sults can be applied to real graphs, finding very good 
agreement of the predicted scaling of the clustering coef- 
ficient C(k) with the square of the maximal eigenvector 
Sk of the matrix Nyr(k : k'),i.e. C(k) ~ S k . Moreover we 
can have a condition for predicting the absence of Hamil- 
tonian cycles for the three-core of Internet and protein- 
protein interaction data. The results indicate that degree 
correlations strongly affect the loop frequency. Further 
study would consider how important are fluctuations of 
the number of loops around this average and would con- 
sider the frequency of other subgraphs in correlated scale- 
free networks. 



VI. ACKNOWLEGMENTS 

The work was supported by EVERGROW, integrated 
project No. 1935 in the complex systems initiative of the 
Future and Emerging Technologies directorate of the 1ST 
Priority, EU Sixth Framework and by EU grant HPRN- 
CT-2002-00319,q. 



Press,Cambridge U.K. ,2004). 
[5] A.-L. Barabasi and R. Albert, Science 286, 509 (1999) 
[6] R. Milo, S. Shen-Orr, S. Itzkovitz, N. Kashtan, D. 

Chklovskii and U. Alon, Science 298, 824 (2002). 
[7] A. Vazquez, R. Dobrin, D. Sergi, J. -P. Eckmann, Z. N. 

Oltvai and A.-L. Barabasi PNAS 101, 17940 (2004). 
[8] S. Janson, T. Luczak and A. Rucinski, Random graphs 



7 



(John Wiley & Sons,New York, 2000). 
[9] G. Bianconi and A. Capocci, Phys. Rev. Lett. 90 ,078701 
(2003). 

[10] G. Bianconi and M. Marsili, J STAT P06005 (2005). 

[11] J. Berg and M. Lassig Phys. Rev. Lett. 89 (2002). 

[12] R. Pastor-Satorras, A. Vazquez and A. Vespignani Phys. 

Rev. Lett.87, 258701 (2001). 
[13] M. E. J. Newman, Phys. Rev. Lett. 89, 208701 (2002). 
[14] E. Ravasz, A. L. Somera, D. A. Mongru, Z. N. Oltvai and 

A.-L. Barabasi, Science 297, 1551 (2002). 
[15] G. Bianconi, G. Caldarelli and A. Capocci, Phys. Rev. E 

71, 066116 (2005). 
[16] S. N. Soffer and A. Vazquez Phys. Rev. E 71, (2005). 
[17] The the Internet datasets are the one collecetd by Uni- 



versity of Oregon Route Views project, NLANR and the 
protein-protein interaction datasets are one listed in DIP 
database. 

[18] E. Marinari, R. Monasson and G. Semerjian, 

cond-mat/0507525 (2005). 
[19] E. Marinari and R. Monasson, J STAT P09004 (2004). 
[20] G Caldarelli, A. Capocci, P. De Los Rios and M. A. 

Muhoz, Phys. Rev. Lett. 89, 258702 (2002). 
[21] M. Boguna and R. Pastor-Satorras, Phys. Rev. E 68 

,036112 (2003). 

[22] K.-L. Goh, B. Kahng and D. Kim Phys. Rev. Lett. 87, 
278701 (2001). 



